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Abstract. We study the thermoelectric transport properties in the three-dimensional Anderson model of 
localization near the metal-insulator transition (MIT). In particular, we investigate the dependence of the 
thermoelectric power S, the thermal conductivity K, and the Lorenz number Lq on temperature T. We first 
calculate the T dependence of the chemical potential fi from the number density n of electrons at the MIT 
using an averaged density of states obtained by diagonalization. Without any additional approximation, 
we determine from /i(r) the behavior of 5", K and Lq at low T as the MIT is approached. We find that a 
and K decrease to zero at the MIT as T ^ and show that S does not diverge. Both S and Lo become 
temperature independent at the MIT and depend only on the critical behavior of the conductivity. 

PACS. 61.43.-j Disordered solids - 71.30.-|-h Metal-insulator transitions & other electronic transitions - 
72.15.Cz Electrical and thermal conduction in amorphous & liquid metals & alloys 



1 Introduction 

The Anderson-type metal-insulator transition (MIT) has 
been the subject of investigation for decades since An- 
derson formulated the problem in 1958 ||l|. He proposed 
that increasing the strength of a random potential in a 
three-dimensional (3D) lattice may cause an "absence of 
diffusion" for the electrons. Today, it is widely accepted 
that near this exclusively-disorder-induced MIT the d. c. 
conductivity a behaves as \E — Ed"^ , where Ec is the crit- 
ical energy or the mobility edge at which the MIT oc- 
curs, and J/ is a universal critical exponent ||^ . Numerical 
studies based on the Anderson Hamiltonian of localiza- 
tion have supported this scenario with much evidence 
^,^,|[|^. In measurements of a near the MIT in semi- 
conductors and amorphous alloys this behavior was also 
observed with varying values of v ranging from 0.5-1.3 [Q, 
^,||. It is currently believed that these different exponents 
are caused by interactions in the system ||l0[] . Indeed, an 
MIT may be induced not only by disorder but also by in- 
teractions such as electron-electron and electron-phonon 
interactions, among others ||ll|. Nevertheless, the exper- 
imental confirmation of the critical behavior of a allows 
the use of the Anderson model in order to describe the 
transition between the insulating and the metallic states 
in disordered systems. 

Besides for the conductivity a, experimental investiga- 
tions can also be done for thermoelectric tran sport p rop- 
erties such as the thermoelectric power S ||,|l^,|l3||, the 
thermal conductivity K and the Lorenz number Lq. The 



behavior of these quantities at low temperature T in dis- 
ordered systems close to the MIT has so far not been sat- 
isfactorily explained. In particular, some authors have ar- 

fd that S diverges or that it remains constant 

[l^ as the MIT is approached from the metallic side. 
In addition, |S'| at the MIT has been predicted [|l^ to 
be of the order of ~ 200/iV/K. On the other hand, mea- 
surements of S close to the MIT conducted on semicon- 
ductors for T < 1 K |13|[_ and on amorphous alloys in the 
range 5 K< T < 350 Kl| yield values of the order of 0.1- 
1 fjY/K. They also showed that S can either be negative 
or positive depending on the donor concentration in semi- 
conductors or the chemical composition of the alloy. The 
large difference between the theoretical and experimental 
values is still not resolved. 

The objective of this paper is to study the behavior of 
the thermoelectric transport properties for the Anderson 
model of localization in disordered systems near the MIT 
at low T. We clarify the above mentioned difference in the 
theoretical calculations for S, by showing that the radius 
of convergence for the Sommerfeld expansion used in Refs. 
@|l|l is zero at the MIT. We show that S* is a finite 
constant at the MIT as argued in Refs. p5ln9]. Besides 
for S, we also compute the T dependence lor a, K, and 
Lq. Our approach is neither restricted to a low- or high-T 
expansion as m Refs. Hill, nor confined to the critical 



regime as in Ref. [16 
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We shall first introduce the model in Sec.|^. Then in 
Sees. ^ and ^ we review the thermoelectric transport prop- 
erties in the framework of linear response and the present 
formulations in calculating them. In Sec. || we shall show 
how to calculate the T dependence of these properties. 
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Fig. 1. The density of states of a 3D Anderson model, aver- 
aged over many disorder realizations with W — 12. The solid 
vertical lines at —Ec and Ec denote the mobility edges. 



Results of these calculations are then presented in Sec.^. 
Lastly, in Sec. we discuss the relevance of our study to 
the experiments. 



2 The Anderson Model of Localization 

The Anderson model pi is described by the Hamiltonian 



(1) 



where is the potential energy at the site j of a regular 
cubic lattice and is assumed to be randomly distributed 
in the range [— T/F/2, WY2] throughout this work. The hop- 
ping parameters tij are restricted to nearest neighbors. For 
this system, at strong enough disorder and in the absence 
of a magnetic field, the one-particle wavefunctions become 
exponentially localized at T = and a vanishes . Illus- 
trating this, we refer to Fig.|^ where we show the density 
of states p(E) obtained by diagonalizingthe Hamiltonian 
(|l|) with the Lanczos method as in Ref. The states 

in the band tails with energy \E\ > Ec are localized within 
finite regions of space in the system at T = When 
the Fermi energy Ep is within these tails at T = the 
system is insulating. Otherwise, if \Ef\ < Ec the system 
is metallic. The critical behavior of a is given by 



a{E) 



1 



E 



0, 



l^^l < Ec, 
\E\ > Ec, 



where ctq is a constant and v is the conductivity exponent 
10. Thus, Ec is called the mobility edge since it sepa- 
rates localized from extended states. At the critical dis- 
order Wc = 16.5, the mobility edge occurs at Ec = 0, 
all states with > are localized |^,^ and states with 
E = are multifractal [^|l^. The value of i/ has been 
computed from the non-linear sigma- model , transfer- 
matrix methods ||^,^, Green functions methods and 
energy-level statistics [^|20|. Here we have chosen v = 1.3, 
which is in agreement with experimental results in Si:P Q 
and the numerical data of Ref. jsj. More recent numeri- 
cal results computed with higher accuracy, suggest 
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Fig. 2. In an open circuit, a temperature gradient VT induces 
an electric field E in the opposite direction which opposes the 
thermal flow of electrons. 



that = 1.5 ± 0.1. As we shall show later, this difference 
only slightly modifies our results. We emphasize that the 
Hamiltonian only incorporates the electronic degrees 
of freedom of a disordered system and further excitations 
such as lattice vibrations are not included. 

For comparison with the experimental results, we mea- 
sure a in Eq. (||) in units of fl~^cm~^. We fix the energy 
scale by setting tij = 1 eV. Hence the band width of Fig. 

tis comparable to the band width of amorphous alloys 
1|. Furthermore, the experimental investigations of the 
thermoelectric power S in arnorphous alloys have been 
done at high electron filling [g2| and thus we will mostly 
concentrate on the MIT at Er- 



3 Linear Thermoelectric Effects 

3.1 Definition of the Transport Properties 

Thermoelectric effects in a system are due mainly to the 
presence of a temperature gradient VT and an electric 
field E We recaU that in the absence of VT with 
E 7^ 0, the electric current density (j) flowing at a point 
in a conductor is directly proportional to E, 



0) = 'tE 



(3) 



By applying a finite gradient VT in an open circuit, elec- 
trons, the thermal conductors, would flow towards the 
low-T end as shown in Fig.|^. This causes a build-up of 
negative charges at the low-T end and a depletion of neg- 
ative charges at the high-T end. Consequently, this sets 
up an electric field E which opposes the thermal flow of 
electrons. For small VT, it is given as 



E = 5VT . 



(4) 



This equation deflnes the thermopower S. In the Som- 
merfeld free electron model of metals, S is found to be 
directly proportional to — T [|2^. Note that the negative 
sign is brought about by the charge of the thermal con- 
ductors. For small VT, the flow of heat in a system is 
proportional to VT. Fourier's Law gives this as 



(j,) - K{-VT) 



(5) 



where (jg) is the heat current density and K is the ther- 
mal conductivity p^ . At low T, the phonon contribution 
to a and K becomes negligible compared to the electronic 
part j23j. As T — + 0, ct approaches a constant and K be- 
comes linear in T. One can then verify the empirical law 
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of Wiedemann and Franz which says that the ratio of K 
and a is directly proportional to T [^|2^. The propor- 
tionality coefficient is known as the Lorenz number Lq, 



(6) 



where e is the electron charge and fcs is the Boltzmann 
constant. For metals, it takes the universal value 7r^/3 p3| , 
|25|| . Strictly speaking, the law of Wiedemann and Franz 
is valid at very low T (< 10 K) and at high (room) T. 
This is because in these regions the electrons are scattered 
elastically. At T ~ 10 — 100 K deviations from the law are 
observed which imply that K/aT depends on T. 

In summary, Eqs. (|^)-(|^) express the phenomenological 
description of the transport properties. 



3.2 The Equations of Linear Response 

A more compact and general way of looking at these ther- 
moelectric "forces" and effects is as follows: the responses 
of a system to E and VT up to linear order Eq] are 



(j) 



(|e|LnE 



and 



(jg) = |er2(|e|L2iE-L22T-iVT) 



(7) 



(8) 



The kinetic coefficients are the keys to calculating the 
transport properties theoretically. Using Ohm's law (0) in 
Eq. (|^), we obtain 

a^Ln- (9) 

Also from Eq. (^, S, measured under the condition of zero 
electric current, is expressed as 



(10) 



\e\TLn 

With the same condition, Eq. (||) yields 
_ L22L11 — L21L12 

Wtl^^ ■ 

Lo is given as 

_ L22L11 — L21L12 
^0 — 



From Eq. 



(11) 



(12) 



Therefore, we will be able to determine the transport prop- 
erties once we know the coefficients Ly . We note that ir 
the absence of a magnetic field, as considered in this work. 



the Onsager relation L21 — L12 holds |26 

Eliminating the kinetic coefRcients in Eqs. (|^) and 
in favor of the transport properties, we obtain 

(j) = crE - aSVT 



and 



T 



SO) 



KVT 
T ' 



(13) 



(14) 



Here, {jq)/T is simply the entropy current density pof . 
Hence, the thermopower is just the entropy transported 
per Coulomb by the flow of thermal conductors. Accord- 
ing to the third law of thermodynamics, the entropy of a 
system and, thus, also {jn)/T willgo to zero as T ^ 0. 
We can check with Eqs. (O) and (nj) that this is satisfied 
by our calculations in the 3D Anderson model. 



3.3 Application to the Anderson Transition 

In general, the linear response coefRcients are obtained 
through the Chester-Thellung-Kubo-Greenwood (CTKG) 
formulation |2^,^. The kinetic coefRcients are expressed 

as 



Ln 



L12 
and 
L22 = 



AiE) 



A{E) [E - ^,{T)] 



dE 



dE. 



(15) 



df{E,fi,T) 



dE 



dE. 



(16) 



A{E) [E-^,{T)f 



df{E,ij,T) 



dE 



dE, (17) 



where A{E) contains all the system-dependent features, 
fi(T) is the chemical potential and 

f{E, T) = 1/ {1 + cxp([i? - ^liT)]/kBT)} (18) 

is the Fermi function. The CTKG approach inherently 
assumes that the electrons are noninteracting and that 
they are scattered elastically by static impurities or by lat- 
tice vibrations. A nice feature of this formulation is that 
all microscopic details of the system such as the depen- 
dence on the strength of the disorder enter only in A{E). 
This function A{E) can be calculated in the context of 
the relaxation-time approximation [p^ . However, an ex- 
act evaluation of is difficult, if not impossible, since it 
relies on the exact knowledge of the energy and T depen- 
dence of the relaxation time. In most instances, these are 
not known. 

In order to incorporate the Anderson model and the 
MIT in the CTKG formulation, a different approach is 
taken: We have seen in Eq. (^ that the d.c. conductivity 
is just Lii. Thus, to take into account the MIT in this 
formulation, we identify A{E) with a{E) given in Eq. (||). 
The L,j in Eqs. (|l5|)-(|l3) can now be easily evaluated close 
to the MIT without any approximation, once the T depen- 
dence of the chemical potential /i is known. Unfortunately, 
this is not known for the experimental systems under con- 
sideration [Q,||,||,|l^,|l3|, nor for the 3D Anderson model. 
Thus one has to resort to approximate estimations of /i, 
as we do next, or to numerical calculations, as we shall do 
in the next sections. 



4 Evaluation of the Transport Coefficients 
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4.1 Sommerfeld expansion in the metallic regime 

Circumventing the computation of /i(T), one can use that 
—df/dE is appreciable only in an energy range of the 
order of ksT near fi Ep- The lowest non-zero T cor- 
rections for the Lij are then accessible by the Sommer- 
feld expansion provided that A{E) is nonsingular and 
slowly varying in this region. Hence, in the limit T — ^ 0, 
the transport properties are Esl 



AiEF) + —ikBTf 





d'^A{E) 



dE^ 



s 



TT^klT dA{E) 



i\e\A{EF) dE 



E — E p 



^ ^ 7:^k%T I ^^^^^ _ n^iksT)' \ dAiE) 



3e2 

and consequently 



3A{Ef) 



dE 



E — Ep 



_ TT^ j n^ksTy [ dA{E) 
° 3 I 'S[A{EfW [ dE 



(19) 



(20) 



(21) 



(22) 



E — Ef 



In the derivations of S, K, and Lq, the term of order 
in Eq. ( p^ ) has been ignored as is customary. We remark 
that the terms of order in Eqs. (|2l|) and ( p^ are usu- 
ally dropped, too. In this case in the metallic regime, Lq 
reduces to the universal value 7r^/3 p^ . 

The above approach was adopted in Refs. |l^ and |^ 
to study thermoelectric transport properties in the metal- 
lic regime close to the MIT. From Eq. (|20|), the authors 
deduce 



S 



lyTT^kgT 

'3\e\{EF-Ec) 



(23) 



In the metallic regime, this linear T dependence of S 
agrees with that of the Sommerfeld model of metals p^ . 
However, setting A{E) = a{E) at the MIT in Eq^ 
is in contradiction to the basic assumption of the Som- 
merfeld expansion, since it is not smoothly varying at 
Ef = Ec. Thus identifying A{E) = a{E) in Eqs. |l| - 12| is 
only valid in the metallic regime with ksT <C \Ec — Ef\- 



4.2 Exact calculation at ^(T) = E^ 

A different approach taken by Enderby and Barnes is to 
fix /i = —Ec at finite T and later take the limit T — > |jl^ . 
Thus, again without knowing the explicit T dependence 
of /i, the coefficients can be evaluated at the MIT. For 
the transport properties they obtain. 



5 = 



aov{kBTyh 
kB ly + l lu+i 



(24) 



(25) 



K 



and 



<JoikBTr+^ 
e^T\EX 



{v + 2)/,+2 



+1 



{v + 2)7^+2 
vl,. 



t\2t2 



(26) 



(27) 



Here h = In 2, = (1 - 2^-'')r{i')C{v) for Re{v) > 
0, V 1, with r{i') and Ci^^) the usual gamma and Rie- 
mann zeta functions. We see that at the MIT, S does not 
diverge nor go to zero but remains a universal constant. Its 
value depends only on the conductivity exponent v. This 
is in contrast to the result (|2^) of the Sommerfeld expan- 
sion. In addition, we find that a oc T'^ and K oc T"^^ 
as T — > 0. Hence, a and K/T approach zero in the same 
way. This signifies that the Wiedemann and Franz law is 
also valid at the MIT recovering an earlier result in Ref. 

obtained via diagrammatic methods. However, at the 
MIT, Lq does not approach 7r^/3 but again depends on i/. 



24)-(|27|) are exact at T values 
16 1 . Thus the T dependence of 
(7, 5, K, and Lq for a given electron density can only be 
determined if one knows the corresponding fi(T). 



Wc emphasize that Eqs 
such that fi{T) - Ec = 



4.3 High-temperature expansion 

In this section, we will study the lowest-order corrections 
to the results obtained before with u(T) = Ec. We do this 
by expanding the Fermi function ( [l8|) for \Ec — fJ-{T)\ ^ 
kBT. In addition, we assume ^i{T) « Ef for the temper- 
ature range considered. This procedure gives 



<JoV{kBTY 

\EX 



En — Ej, 



knT 



(28) 



For the thermopower, the leading-order correction can be 
obtained without expanding f{E,iJ,,T) in Lu and L12. 
This yields a constant for S at the MIT [O . We obtain 



_ kB 



II u 



+1 



Ec — Ef 
kBT 



(29) 



For K and Lq, we again have to use the expansion of 
f{E,fi,T) as in (28) in order to get non-trivial terms. 
The resulting expressions are cumbersome and we thus 
refrain from showing them here. We remark that the basic 
ingredients used in the high-T expansion are somewhat 
contradictory, namely, the expansion is valid for high T 
such that \Ec — Ef\ <C fc^T, whereas /^(T) = Ef is true 
only for T = 0. 

At present, wc thus have various methods of circum- 
venting the explicit computation of fJ-(T). However, their 
ranges of validity are not overlapping and it is a priori not 
clear whether the assumptions for /i(T) are justified for S 
or any of the other transport properties close to the MIT. 
In order to clarify the situation, we numerically compute 
//(T) in the next section and then use the CTKG for- 
mulation to compute the thermal properties without any 
approximation. 
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5 The Numerical Method 

In Eqs. ([l5|)-([l7|), the explicit T dependence of the coeffi- 
cients Lij occurs in f{E,fi,T) and /i(r). More precisely, 
knowing fJ-(T), it is straightforward to evaluate the Lij. 
We recall that, for any set of noninteracting particles, the 
number density of particles n can be determined as 



dEp{E)f(E,ix,T) 



(30) 



where p{E) is again the density of energy levels (in the 
unit volume) as in Fig , [l l. Vice versa, if we know n and 
p{E) we can solve Eq. (|3q ) for p.{T). The density of states 
p{E) for the 3D Anderson model has been obtained for 
different disorder strengths W as outlined in Sec. ^. We 
determine p{E) with an energy resolution of at least 0.1 
meV (~ 1 K). Using p{E), we first numerically calculate 
n at T = for the metallic, critical and insulating regimes 
using the respective Fermi energies \Ef\ < Ec, Ep — Ec, 
and \Ef\ > Ec- With ^ = Ep, we have 



n{Ei, 



dEp{E) 



(31) 



Next, keeping n fixed at n{Ep), we numerically determine 
/z(T) for small T > such that IniEp) — n{p, T)\ is zero. 
Then we increase T and record the respective changes in 
/x(T). Using this result in Eqs. (|l5l)-(|l^) in the CTKG for- 
mulation, we compute Lij by numerical integration and 
subsequently determine the T dependent transport prop- 
erties %-^). 

We consider the disorders W = 8, 12, and 14 where 
we do not have large fluctuations in the density of states. 
These values are not too close to the critical disorder Wc, 
so that we could clearly observe the MIT of Eq. (H) . The 
respective values of E^ have been calculated previously Q 
to be close to 7.0, 7.5, and 8.0. Within our approach, we 
choose Ec to be equal to these values. 



6 Results and Discussions 

Here we show the results obtained for W = 12 with Ec — 
7.5. The results for a, K, and Lq are the same at —Ec and 
Ec since they are functions of in, L22 and L12, only. On 
the other hand, this is not true for S. 



6.1 The Chemical Potential 

In Fig.^, we show how p{T) behaves for the 3D Anderson 
model at Ep—Ec — 0, and ±0.01. To compare results from 
different energy regions we plot the difference of /^(T) from 
Ep. We find that p{T) behaves similarly in the metallic 
and insulating regions and at the MIT for both mobility 
edges at low T. In all cases we observe /i(T) (x . Fur- 
thermore, we see that p{T) at —Ec equals —p{T) at Ec. 
This symmetric behavior with respect to Ep — /i reflects 
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Fig. 3. The temperature dependence of the chemical poten- 
tial /i measured with respect to the Fermi energy near both 
mobility edges. Also shown is p{T) for a free electron gas. The 
solid line denotes p{T) of Eq. 



the symmetry of the density of states at i? = as shown 
inFig.|l|. 

For comparison and as a check to our numerics, we 
also compute with our method p.{T) of a free electron gas. 
The density of states is 



p{E) 



3 n 



2Ep\E 



E 



1/2 



(32) 



and we again use Ep — Ec — 7.5. We remark that this 
value of the mobility edge is in a region where p{E) in- 
creases with E in an analogous way as p{E) for the An- 
derson model at —Ec ■ Thus, as shown in Fig. ^, p{T) of a 
free electron gas is concave upwards as in the case of the 
Anderson model at —Ec. We also plot the result for piT) 
obtained by the usual Sommerfeld expansion for Eq. (pO|), 



Ep - p(T) 



Ep ( irkBT 
2Ep 



(33) 



We see that our numerical approach is in perfect agree- 
ment with the free electron result. 



6.2 The d.c. Conductivity 

In Fig. H we show the T dependence of a. The values of 
Ep we consider and the corresponding fillings n are given 
in Tab. ^ The conductivity at T = remains finite in 
the metallic regime with ct/ctq = |1 — Ep/Ec\'^ , because 
i-df/dE) 5{E - Ep) in Eq. (|l|) as T ^ 0. Cor- 
respondingly, we find a — Q in the insulating regime at 
T = 0. In the critical regime, (t(T ^ 0) ^ y, as derived 
in Rcf. |16 , sec Eq. (p^). We note that as one moves away 
from the critical regime towards the metallic regime one 
finds within the accuracy of our data that a ^ T^. We ob- 
serve that in the metallic regime a increases for increasing 
T. This is different from the behavior in a real metal where 
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tr decreases with increasing T. However, as explained in 
Sec. H, the behavior of a in Fig.jJ is due to the absence of 
phonons in the present model. 

We also show in Fig.^ results of the Sommerfeld ex- 
pansion ( 19) and the high-T expansion ( p8| ) for a. Paradig- 
matic for what is to follow we see that the radius of conver- 
gence of the Sommerfeld expansion decreases for Ep Ec 
and in fact is zero in the critical regime. On the other 
hand, the high-T expansion is very good in the critical 
regime down to T = at i?c = Ep. The small systematic 
differences between our numerical results and the high-T 
expansion for large T are due to the differences in /i(T) 
and Ep. The expansion becomes worse both in the metal- 
lic and insulating regimes for larger T. All of this is in 
complete agreement with the discussion of the expansions 
in Sec. |. 



6.3 The Thermopower 

In Fig. H, we show the behavior of the thermopower at low 
T near the MIT. In the metallic regime, we find 5* ^ 
as T — > 0. At very low T, S oc T as predicted by the 
Sommerfeld expansion (^). We see that the Sommerfeld 
expansion is valid for not too large values of T. But upon 
approaching the critical regime, the expansion becomes 
unre liab le similar to the case of the d.c. conductivity of 
Sec. |6.2| . This behavior persists even if we include higher 
order terms in the derivation of S such as the O(T^) term 
of Eq. (19) as shown in Fig.||. 

Before discussing the critical regime in detail, let us 
turn our attention to the insulating regime. Here, S be- 
comes very large as T 0. We have observed that it 
even appears to approach infinity. A seemingly divergent 
behavior in the insulating regime has also been observed 
for Si:P where it has been attributed to the ther- 
mal activation of charge carriers from Ep to the mobility 
edge Ec- However, there is a simpler way of looking at this 
phenomenon. We refer again to the open circuit in Fig.^. 
Suppose we adjust T at the cooler end such that VT re- 
mains constant. As T ^ both a and K vanish in the case 
of insulators — for K we show this in the next section. 



T (meV) 

4 6 8 10 



Table 1. Differences of Ep and n{Ep) with respect to the 
mobility edge at Ec = 7.5. The density at Ec corresponds to 
n = 97.768%. 
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Fig. 4. The low temperature behavior of the d.c. conductiv- 
ity (J. The symbols are as shown in Tab. |^. The dashed lines 
represent the Sommerfeld expansion result for (j{T) as given 
in Eq. (|is|). For all 8 choices of Ep — Ec, the corresponding 
high-T expansion torn is indicated by solid lines. 



This implies that as T decreases it becomes increasingly 
difficult to move a charge from T to T -I- 5T. We would 
need to exert a larger amount of force, and hence, a larger 
E to do the job. From Eq. (P, this implies a larger S value. 

In the critical regime, i.e., setting Ep = E^ we observe 
in Fig. P that for T ^ the thermopower S approaches a 
value of 228.4 iiY/K. This is exactly the magnitude pre- 
dicted [[6) by Eq. (||) for = 1.3. In the inset of Fig. |, we 
show that the T dependence of S is linear. The nondiver- 
gent behavior of S clearly separates the metallic from the 
insulating regime. Furthermore, just as for cr, the Som- 
merfeld expansion for S breaks down at Ep = Ec, i.e., 
the radius of convergence is zero. Thus, the divergence of 
Eq. ( p^ ) at Ep = Ec reflects this breakdown and is not 
physically relevant. On the other hand, the high-T ex- 
pansion 1 15 nicely reflects the behavior of S close to the 
critical regime as also shown in Fig. ||. For Ep — Ec, the 
high-T expansion (|2^) assumes a constant value of S for 
all T due to setting /x(T) = Ep. This is approximately 
valid, the differences are fairly small as shown in the inset 
of Fig.|. 

We stress that there is no contradiction that S > 
in our calculations whereas 5 < in Ref. ||l^. In Fig.^, 
we compare 5* in energy regions close to Ec and to ~Ec 
Clearly, they have the same magnitude but 5 < 
at —Ec and S > at Ec- The two cases mainly differ in 
their number density n. At —Ec the system is at low filling 
with n = 2.26% while at Ec the system is at high filling 
with n = 97.74%. The sign of S implies that at low filling 
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Fig. 5. The low temperature behavior of the thermopower S. 
The symbols are as shown in Tab. |l[ The dashed lines represent 
the behavior of 5(7) in the metallic regime as given in Eq. (ES) . 
The dot-dashed lines indicate S, calculated with the 0(T^) 
term of Eq. (|l|), for Ep - Ec = -0.01 eV (o) and -0.001 eV 
(o). Solid lines are obtained from the high-T expansion (p9|). 
The inset shows the behavior at Ep — E^ on an enlarged scale. 
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Fig. 6. An example that the magnitude of S{T) is the same in 
metallic regions close to —Ec (■) and Ec (o). The +- symbols 
indicate \S\ for —Ec and \Ep — EJ — 0.01 eV in all cases. 



the thermoelectric conduction is due to electrons and we 
obtain the usual picture as in Fig.H where the induced 
field E is in the direction opposite to that of VT. At high 
filling, S > means that E is directed parallel to VT. This 
can be interpreted as a change in charge transport from 
electrons to holes. We remark that this sign reversal also 
occurs in the insulating as well as in the critical regime. 
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I. 16 K = 0.1 meV 
4.34 K = 0.4 meV 
8.12 K = 0.7 meV 

II. 6 K= 1.0 meV 
22.3 K = 2.0 meV 
S = 228.4 |xV/K 
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Fig. 7. The data of S in Fig. ^ shown as a function of fj, mea- 
sured from Ec — 7.5 eV. The horizontal line indicates the fixed 
point MIT value as given in Eq. (^^. The thin dashed lines 
represent isotherms of 5* calculated using the same method as 
in Rcf. jl6[. The solid line is an isotherm of S obtained from 
Eq. (ph for T = 22.3 K. 
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Fig. 8. Scaling plot of the thermopower S. The thick dashed 
line indicates the fixed point value at the MIT, the solid line 
represents the high-T expansion (|2^) , and the thin dashed line 
shows the Sommerfeld expansion. The inset shows the differ- 
ence in the scaling when plotting S for Ep — Ec = —0.001 eV 
as function of Ec) /ksT (open symbols) or [Ep — Ec) /ksT 
(filled symbols). 



In Fig. in we take the data of Fig.|| and plot them as a 
function of /i — i?c. Our data coincides with the isothermal 
lines which were calculated according to Ref. |l6[ by nu- 
merically integrating L12 and Ln for a particular T to get 
S. We observe that all isotherms of the insulating (/x > E^) 
and the metallic (/it < Ec) regimes cross at n = Ec and 
5 = 228.4 ^V/K. Comparing with Eq. (jH), we again find 
that the Sommerfeld expansion does not give the correct 
behavior of S in the critical regime. 

The data presented in Fig.|^ suggest that one can scale 
them onto a single scaling curve. In Fig.||, we show that 
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this is indeed true, when plotting S* as a function of (fi — 
Ec)/kBT. We emphasize that the scaling is very good and 
the small width of the scaling curve is only due to the 
size of the symbols. The result for the high-T expansion 
is indicated in Fig.|| by a solid line. It is good close to 
the MIT. In the metallic regime, the Sommerfeld expan- 
sion correctly captures the decrease of S for large negative 
values of (/i — Ec)/kBT. We remark that a scaling with 
{Ep — Ec)/kBT as predicted in Ref. is approximately 
valid. The differences are very small as shown in the inset 
of Fig.H. 



6.4 The Thermal Conductivity and the Lorenz Number 

In Fig. ^ we show the T dependence of the thermal con- 
ductivity K. We see that if — > as T ^ whether it be 
in the metallic or insulating regime. We note again that 
this simple behavior is due to the fact that our model does 
not incorporate phonon contributions. The T dependence 
of K varies whether one is in the metallic regime or in 
the insulating regime and how far one is from the MIT. 
Directly at the MIT, we find that X ^ as T^+i confirm- 
ing the T dependence of K as given in Eq. (p6|). Near the 
localization MIT, the T dependence of K/T is thus the 
same as for a in agreement with Ref. Again, we see 
that the Sommerfeld expansion ( pi] ) is reasonable only at 
low T in the metallic regime. As for a and S, we see that 
the high-T expansion is again fairly good in the vicinity 
of the critical regime. 

At this point we are able to determine the behavior 
of the entropy in the system as T — > 0. In the metallic 
regime, S and K vanish as T — > 0, while in the critical 
and insulating regime, a and K vanish as T ^ 0. Ap- 
plying these results to Eqs. ( |l3|) and (^4|) yields that for 
all regimes the entropy current density (jq) /T vanishes as 
T ^ 0. Therefore, we find that the third law of thermo- 
dynamics is satisfied for our numerical results of the 3D 
Anderson model. 

Next, we present the Lorenz number (|^) as a func- 
tion of T in Fig. no. In the metallic regime, we obtain the 
universal value 7p/3 as T 0. Note that for a metal 
this value should hold up to room T However, our 
results for the Anderson model show a nontrivial T de- 
pendence. One might have hoped that the higher-order 
terms in Eq. ( |2^ ) could adequately reflect the T depen- 
dence of our Ln data. However, this is not the case as 
shown in Fig. nG. This indicates that even if we incorpo- 
rate higher order T corrections the Sommerfeld expansion 
will not give the right behavior of Lq near the MIT. We 
emphasize that the radius of convergence of Eq. (|2^) is 
even smaller than for a, S and K. Similarly, the high-T 
expansion is also much worse than previously for a, S and 
K. Thus in addition to the results for the critical regime, 
we only show in Fig.jl^ the results for nearby data sets 
in the insulating and metallic regimes. The T dependence 
of Lq is linear as shown in the inset of Fig. ^ As before 
for S, the high-T expansion does not reproduce this. At 
the MIT, Lq = 2.4142. This is again the predicted ^ 
i^-dependent value as given in Eq. (p7|). 
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Fig. 9. The thermal conductivity K as a function of temper- 
ature. The symbols are as shown in Tab. The dashed lines 
were obtained in 0(r) from the Sommerfeld expansion (|2l| ) 
for the metallic regime. The results of the high-T expansion 
for the 8 choices of Ep ~ Ec are indicated by solid lines. 
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Fig. 10. The Lorenz number Lo as function of temperature. 
The symbols are as shown in Tab. ^. The dashed circles mark 
the values of Lq at T = for metallic and insulating regimes. 
The dashed lines were obtained from Eq. (^). The results of 
the high-T expansion for Ef-Ec = eV, ±0.001 eV and 0.003 
eV are indicated by solid lines. The inset shows the behavior 
at Ep — Ec on an enlarged scale. 
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In the insulating regime, one can show analytically by 
taking the appropriate limits that Lq approaches + 1 as 
T ^ 0. In agreement with this, we find that Lq — 2.3 at 
T = in Fig. At first glance, it may appear surprising 
that a transport property in the insulating regime could be 
determined by a universal constant of the critical regime 
such as v. However, in the evaluation of the coefficients 
Lij^ the derivative of the Fermi function for any finite T 
decays exponentially and thus one will always have a non- 
zero overlap with the critical regime. In the evaluation of 
Eq. (p^), this v dependence survives in the limit T — > 0. In 
real materials, we expect the relevant high-energy transfer 
processes to be dominated by other scattering events and 
thus Lq should be different. Nevertheless, for the present 
model, this v dependence holds. 



6.5 Possible Scenarios in the Critical Regime 

The results presented in Sec. |6.3| for the thermopower at 
the MIT show that S = 228.4 ^V/K for = 1.3. This 
value is 2 orders of magnitude larger than those measured 
near the MIT |,|l2|, 
troduction, the cone 
periments are either 
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However, as mentioned in the in- 
uctivity exponents found in many ex- 
close v — 0.5 or to 1 and one 
might hope that this difference may explain the small ex- 
perimental value of S. Also, recent numerical studies of 
the MIT by transfer-matrix methods together with non- 
linear finite-size scaling find — 1.57 ± 0.03 In Tab. 
^ we summarize the values of S and Lq at the MIT for 
these conductivity exponents. We see that all S values still 
differ by 2 orders of magnitude from the experimental re- 
sults. Furthermore, we note that our results for S and Lq 
are independent of the unit of energy. Even if, instead of 
1 eV, we had used ttj = 1 me V, which is appropriate in 
the doped semiconductors M|9[|l^,^, we would still ob- 
tain the values as in Tab. 2f Thus our numerical results 
for the thermopower of the Anderson model at the MIT 
show a large discrepancy from experimental results. This 
may be due to our assumption of the validity of Eq. (^) for 
a large range of energies, or due to the absence of a true 
Anderson-type MIT in real materials, or due to problems 
in the experiments. 

A different scenario for a disorder driven MIT has been 
proposed by Mott, who argued that the MIT from the 
metallic state to the insulating state is discontinuous p3 . 



Table 2. The thermopower and the Lorenz number at the 
MIT for a 3D Anderson model evaluated for various f at Ec = 
7.5 eV. The values for v — 0.5 and 1 have already been shown 
in Ref. M. 



V 


S 






(mV/K) 




0.5 


163.5 


1.7761 


1.0 


204.5 


2.1721 


1.3 


228.4 


2.4142 


1.57 


249.7 


2.6372 



Results supporting such a behavior have been found ex- 
perimentally jri|,p3[. According to this scenario, a drops 
from a finite value cr„ij„ to zero for T = at the MIT. 
This minimum metallic conductivity Omin was estimated 
by Mott to be 

1 

Omm — - ^ (34) 
a n 

where a is some microscopic length of the system such as 



the inverse of the Fermi wave number, a 



As sum- 



marized in Ref. |11| , experiments in non-crystalline mate- 
rials seem to indicate that amin > 300 f2~^cm~^. Let us 
assume the behavior of <j{E) close to the MIT to be 



a{E) = 



0, 



1^1 < 
\E\ > Er, 



(35) 



with a„iin — 300 17 ^cm ^. Using the numerical approach 
of Sec. |, we obtain S = 119.5 /iV/K at the MIT. This 
value is still rather large and thus the assumption of a 
minimum metallic conductivity as in Eq. ( p5| ) cannot ex- 
plain the discrepancy from the experimental results. We 
remark that the order of magnitude of S is not changed 
appreciably, even if we add to the metallic side of Eq. (|3^ ) 
a term as given in Eq. with ctq a- few hundred Q~^cm~^ 
and v = 1. 

Lastly, we note that the transport properties calcu- 
lated for = 8 and 14 do not differ from those obtained 
for W — 12 in both the metallic and insulating regions 
provided we are at temperatures T < lOOK. For S and 
Lq at the MIT we obtain the same values as for W = 12. 
Again we observe that both S and Lq approach these val- 
ues linearly with T, but with different slopes. Our results 
show that the higher the disorder strength the smaller the 
magnitude of the slope. 



7 Conclusions 

In this paper, we investigated the thermoelectric effects 
in the 3D Anderson model near the MIT. The T depen- 
dence of the transport properties is determined by fi{T). 
We were able to compute /x(T) by numerically invert- 
ing the formula for the number density n(fi, T) of non- 
interacting particles. Using the result for /^(T), we calcu- 
lated the thermoelectric transport properties within the 
Chester-Thellung-Kubo-Greenwood formulation of linear 
response. As T — > in the metallic regime we verified that 
a remains finite, 5 — > 0, -fC ^ and Lq 7r^/3. On the 
other hand, in the insulating regime, S — > oo. This we at- 
tribute to both a and K going to zero. Thus, it becomes 
increasingly difficult to achieve equilibrium and, hence, 
the system requires E ^ oo. For Lq, we obtained a uni- 
versal value of 1 even in the insulating regime. Directly 
at the MIT, the thermoelectric transport properties agree 
with those obtained in Ref. [|l6|. Namely, as T — > 0, we 
found <7 '^T", K ^ T-^+i, while Lq const. 

The thermopower S also remains nearly constant in 
the critical regime and, in particular, it does not diverge 
at the MIT in contrast to earlier calculations using the 
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Sommerfeld expansion at low T [Q . Here we showed that 
the difference is not so much due to an order of Umits 
problem, but rather reflects the breakdown of convergence 
of the Sommerfeld expansion at the MIT . Our result 
is supported by scaling data for S at different values of 
T and Ep onto a single curve which is continuous across 
the transition. Some of the experiments for S have 

been influenced by the Sommerfeld expansion such that 
the authors plot their results as S/T . We remark that in 
such a plot the signature of the MIT is hard to identify, 
since S/T at the MIT diverges as T ^ solely due to the 
decrease in T. Our results suggest that plots as in Figs. || 
and @ should show the MIT more clearly. 

The value of S is at least two orders of magnitude 
larger than observed in experiments |^,^|l^. This large 
discrepancy may be due to the ingredients of our study, 
namely, we assumed that a simple power-law behavior of 
the conductivity (j{E) as in Eq. (g) was valid even for 
E <^ Ec and E ^ Ec- Furthermore, we assumed that it 
is enough to consider an averaged density of states p{E). 
While the first assumption is of course crucial, the sec- 
ond assumption is of less importance as we have checked: 
Local fluctuations in p{E) will lead to fluctuations in the 
thermoelectric properties for flnite T, but do not lead to a 
different T — ^ behavior: S remains finite with values as 
given in Tab. |^. Moreover, averaging over many samples 
yields a suppression of these fluctuations and a recovery 
of the previous behavior for finite T. In this context, we 
remark that — naively assuming all other parts of the 
derivation are unchanged — implications of many-particle 
interactions such as a reduced single-particle density of 
states at Ep [Q, will only modify the T dependence of 
/K. Consequently, the T dependencies of S*, a, K, and Lq 
may be different, but their values at the MIT remain the 
same. 

Our results also suggest that the critical regime is very 
small. Namely, as the filling increases slightly from n — 
97.74% to 97.80%, the behavior of the system changes 
from metallic to critical and finally to insulating. Up to the 
best of our knowledge, such small changes in the electron 
concentration have not been used in the measurements of 
S as in Refs. |||,|l^,^. We emphasize that such a fine 
tuning of n is not essential for measurements of a as is 
apparent from Fig. ^. 

Of course, one may also speculate that these re- 
sults suggest that a true Anderson-type MIT has not yet 
been observed in the experiments. 
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